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ABSTRACT 

The vast majority of dwarf satellites orbiting the Milky Way and M31 are quenched, 
while comparable galaxies in the field are gas-rich and star-forming. Assuming that this 
dichotomy is driven by environmental quenching, we use the ELVIS suite of A-body 
simulations to constrain the characteristic timescale upon which satellites must quench 
following infall into the virial volumes of their hosts. The high satellite quenched 
fraction observed in the Local Group demands an extremely short quenching timescale 
2 Gyr) for dwarf satellites in the mass range M* ~ 10® — 10® M©. This quenching 
timescale is significantly shorter than that required to explain the quenched fraction of 
more massive satellites (~ 8 Gyr), both in the Local Group and in more massive host 
halos, suggesting a dramatic change in the dominant satellite quenching mechanism 
at M* < 10® Mq. Combining our work with the results of complementary analyses in 
the literature, we conclude that the suppression of star formation in massive satellites 
(M* ~ 10® — 10^^ Mq) is broadly consistent with being driven by starvation, such that 
the satellite quenching timescale corresponds to the cold gas depletion time. Below 
a critical stellar mass scale of ~ 10® Mq, however, the required quenching times are 
much shorter than the expected cold gas depletion times. Instead, quenching must 
act on a timescale comparable to the dynamical time of the host halo. We posit that 
ram-pressure stripping can naturally explain this behavior, with the critical mass (of 
M* ^ 10® Mq) corresponding to halos with gravitational restoring forces that are too 
weak to overcome the drag force encountered when moving through an extended, hot 
circumgalactic medium. 
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- galaxies: star formation 


1 INTRODUCTION 


Foremost among the results of galaxy surveys over the last 
decade has been the realization that the galaxy population 
at 2 < 2 is bimodal in nature (e.g. |Strateva et al.||2001 


Baldry et al. 2004 Bell et al. 2004 


Cooper et al. 20061. 


That is, galaxies both locally and out to intermediate red- 
shift are effectively described as one of two distinct types: 
red, early-type galaxies lacking significant star formation 
and blue, late-type galaxies with active star formation. In 
color-magnitude space, the red galaxies populate a tight re¬ 
lation (often called the red sequence), while the distribu¬ 
tion of blue galaxies is more scattered (sometimes referred 
to as the blue cloud). While the red and blue populations 
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comprise approximately equal portions of the cosmic stellar 
mass budget at 2 ~ 1, galaxies on the red sequence domi¬ 
nate today, following a growth in stellar mass within the red 
population of roughly a factor of 2 over the past 7 Gyr (Bell] 


et al. 12004] |Bundy et al.|2006| [Faber et al.|2007| [Brown et*^ 


20071. Despite uncertainty regarding the particular physical 
process(es) at play, the suppression (or quenching) of star 
formation in blue galaxies, thereby making them red, is one 
of the principal drivers of this dramatic growth in the num¬ 
ber density of quiescent systems at late cosmic time. 


At both low and intermediate redshift, the local envi¬ 
ronment of a galaxy is known to be well-correlated with the 
suppression of star formation, such that passive or quies¬ 
cent galaxies preferentially live in higher-density environ¬ 


ments i Balogh et al.||2004 Kauffmann et al~]|2004 Blanton] 
et al.|2005 Cooper et al.|2006 2007| 2010aI. While the stel- 
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2 Fillingham et al. 


lar mass or surface mass density of a galaxy may be more 


closely connected to quenching for massive systems (Peng 
|et al.|[2010[ [Cheung et al.|[2012[ [Woo et al.||2013| ), recent 
work has shown that environment is likely the dominant 
driver of quenching at the lowest mass scales (M* < 10® Mq , 
[Geha et al.||201!^ . For example, studies comparing satellite 
galaxies to isolated field systems of similar stellar mass in 
the local Universe find that satellites tend to exhibit lower 
star formation rates, more bulge-dominated morphologies, 
as well as older and more metal-rich stellar populations 


([Baldry et al.|2006 

van den Bosch et al.|2008 

Cooper et al. 

2010bt jPasquali et al.||20101 [Tolferud et al.| 

20111 [Phillips 


satellite galaxies is commonly referred to as “environmen¬ 
tal quenching”. At present, galaxy formation models fail to 
reproduce this environment-dependent effect in detail, with 
both semi-analytic and hydrodynamical models overpredict- 


ing the number of quenched satellites at low masses ( 

Kimm 

et al.j 

2009||Weinmann et al.||2011| |2012||Hirschmann et al. 

2014 

Wang et al.|2014 

1 . 


Understanding the physics of environmental quenching 
is complicated by the wide range of physical processes that 
could be responsible for suppressing star formation in low- 
mass satellites. For example, various quenching mechanisms 
preferentially operate on satellite galaxies within overdense 
environments, including [i] “starvation” or “strangulation” - 
whereby gas accretion onto a satellite galaxy is halted follow¬ 
ing infall, thus eventually eliminating the fuel for star forma¬ 
tion ([Larson, Tinsley fc Caldwell||198d||Kawata fc Mulchaeyj 


20081, [ii\ “harassment” - by which close encounters be¬ 
tween densely packed cluster or group members strip gas 
from around the interacting galaxies ( [Moore et al.[[T996[ ), 
and [in] “ram-pressure stripping” - where the cold dense 
gas at the center of a satellite (or the hot halo surrounding 
it) is removed from the galaxy as a result of a high-speed 


interaction with the hot gas halo of the host (Gunn & Gott 
lOT^lBekkipOO^ . 


To further complicate matters, each of these processes 
may be important in different host halos and at different 
satellite mass scales. In an attempt to differentiate between 
these various physical mechanisms, we aim to measure the 
timescale upon which satellite quenching occurs across a 
broad range of satellite stellar mass. Throughout this work, 
we define this quenching timescale (rquench) relative to the 
time at which a satellite was accreted onto its current host 
system (i.e. relative to the infall time of the satellite). While 
quenching via starvation will proceed according to the avail¬ 
able cold gas reservoir for star formation within an ac¬ 
creted satellite, ram-pressure stripping will largely act on 
a timescale set by the density distribution of the host’s hot 
halo and infall velocity of the satellite. As such, constraining 
the timescale for satellite quenching may therefore serve as a 
critical step towards differentiating between various physical 
mechanisms. 

By comparing satellite samples at low and intermedi¬ 
ate redshift to galaxy and dark matter-only simulations, 
several recent analyses of satellites at high stellar mass 
(> 5 X 10® Mq) point towards relatively long quenching 
timescales of several Gyr or more (Wetzel et al.|2013 De Lu- 


cia et al.|2012 1. Remarkably, whe n pushing this analys is to 
lower masses (~ 10®'® —10®'® Mq), Wheeler et al. (20141 find 


that quenching remains surprisingly inefficient, such that the 


average quenching timescale increases to > 7 Gyr. Here, we 
extend these previous analyses of quenching timescales to 
yet lower masses by studying the satellite galaxies of the 
Local Group. Through comparison to simulations, we aim 
to establish a coherent picture of the timescale upon which 
satellite quenching acts, spanning more than five orders of 
magnitude in satellite mass. In we detail the popula¬ 
tion of Local Group satellites along with an estimation of 
the observed quenched fraction. In m and 0 we describe 
the relevant A^-body simulations and physically-motivated 
quenching models, respectively, detailing the relevant con¬ 
straints on the derived quenching timescale at low masses. 
Finally, in i|^and ^ we discuss the implications of our anal¬ 
ysis, including a “back of the envelope” explanation of the 
underlying physics operating on low-mass satellites and a 
summary of our main results. 


2 LOCAL GROUP SATELLITES 


In assembling our sample of Milky Way and M31 satellites, 
we utilize the compilation of local dwarf galaxy properties 
presented by [McConnachie[ ( [2012[ ) . To minimize the poten¬ 
tial impact of incompleteness at low masses, we limit our 
analysis to those dwarfs with measured stellar masses in the 
range 10® — 10® Mq; however, as discussed in more detail 
below, including systems at lower masses has very little im¬ 
pact on our qualitative or quantitative results. From this 
population of local dwarfs, satellite systems are selected to 
be within 300 kpc of the Milky Way or M31 and bound 
with respect to their host — i.e. %/3 • Vjos < Uesc(r), where 
Vesc{r) is the escape velocity at a radial distance of r for an 
NFW halo with a concentration (c) of 8 and a total halo 
mass of 2 X 10^® To account for potential tangen¬ 

tial motion, we mul t iply a ll observed line-of-sight velocities 


from 


McGonnachie 


( |2012[ ) by a factor of \/3 in evaluating 
whether each system is bound. Within the stellar mass range 
of 10® Mq to 10® Mq, Figureshows the position and ve¬ 
locity of all Local Group dwarfs relative to their host (or 
nearest massive neighbor in the case of field systems such as 
the Pegasus dirr). 

We identify 12 satellite systems meeting our selection 
criteria (see Table 1), with the final set of satellites and cor¬ 
responding quenched fraction relatively independent of the 
adopted selection criteria. For instance, field dwarfs such as 
NGC 6822, Leo T, and Phoenix, which are nearly in our 
stellar mass limits, are located significantly beyond our ra¬ 
dial selection limit (i.e. > 400 kpc from the Milky Way). 
Moreover, while the observed line-of-sight velocity of Leo I 
suggests that it may be unbound (see Fig.[^, recent proper 
motion observations conclude that the system is very likely 
bound to the Milky Way even at a distance of more than 
250 kpc ( Sohn et al.[[2013 Boylan-Kolchin et al. [[2013 1. Fi¬ 
nally, given the high metallicity and lack of overdensity in 
RR Lyrae stars, we categorize Canis Major as an overden¬ 
sity in the Milky Way and not a satellite (Momany et al. 
2004| [Butler et al.||2007[ [Mateu et al.||2009[ buF see also 
Martfnez-Delgado et al.|2005|. 


^ Varying the adopted host halo mass and concentration (and 
thus our satellite selection limits) does not have a significant effect 
on our final sample of Local Group satellites. 
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Figure 1. Host-centric velocity versus distance for all known Lo¬ 
cal Group galaxies in the stellar mass range of 10® — 10® Mq. To 
account for unknown tangential velocities, the observed line-of- 
sight velocities have been multiplied by a factor of VS- The dotted 
vertical line at 300 kpc corresponds to the expected virial radius 
of the Milky Way and M31 and represents the radial limit for 
our satellite selection. The dashed lines show the escape velocity 
for an NFW halo with a concentration (c) of 8 and a total mass 
of 2 X 10^^ Satellites are color-coded as star-forming (blue 

circles) or quenched (red diamonds) according to their observed 
Hi gas fraction (see Fig.[^. 


To determine whether a given satellite is quenched, we 
use the observed Hi gas fraction as a proxy for the current 
star formation rate (SFR), defining dwarfs with atomic gas 
fractions less than 10% to be quenched (see Fig. [^. This 
criterion leads to only one satellite in our sample, IC 10, 
categorized as currently star-forming and yields a quenched 
fraction of > 90% for satellites bound to the Milky Way 
and M31 in the selected mass range. Alternative definitions 
of quenched versus star-forming - e.g. using star-formation 
histories derived from photometric observations of spatially- 


Figure 2. Observed Hi gas fraction as a function of distance 
from the nearest host galaxy for local dwarfs included in our mass 
range of 10® — 10® Mq . Systems with an atomic gas fraction less 
than 10% are classified as quenched (red diamonds), with more 
gas-rich dwarfs identified as star-forming (blue circles). At these 
low masses, roughly 90% of satellites are quenched in the Local 
Group. 


sion with the disk of both the Milky Way and M31, the 
focus of our analysis is the relative fraction of star-forming 
versus quenched systems. Assuming that neither quenched 
nor star-forming systems preferentially reside in or behind 
the disk, the quenched fraction should be independent of 
any such incompleteness. From an analysis of the Palomar 
and UK Schmidt photographic plates, which cover more 
than 20,000 square degrees of sky at high Galactic latitude 
(|6| > 20°), |Irwin| ( |1994| ) finds that the known satellite pop¬ 
ulation for the Milky Way is complete down to below our 
adopted mass limit 3 x 10® Mq). This conclusion is 
supported by deeper imaging data from the Sloan Digital 


resolved stellar populations (Weisz et al. 2014a 20151 or 

et al. 20101, and Dark Energy Survey (DES, Diehl et al. 

star-formation rates inferred from Ha emission (Kaisin & 

2014P , which cover a significant portion of the northern and 


southern sky and are complete to beyond the host virial ra- 


similar results. Moreover, our measurement of the satellite 
quenched fraction is not strongly dependent on the adopted 
mass range, as the quenched fraction is largely unchanged 
if we extend the lower mass limit down to 3 x 10® Mq 
so as to account for potential uncertainties or scatter in 
the stellar mass measurements. Over this expanded stel¬ 
lar mass range, the sample would then also include Leo II, 
And XXI, And XXV, and LGS 3. While Hi measurements 
suggest that LGS 3 is likely star-forming, a variety of ob¬ 
servations, including star formation rates inferred from Ha 
narrow-band imaging, indicate that Leo H, And XXI, and 


2008| [Koposov et al.||2008| [Laevens et al.||2014| [Koposov 

et al.|2015 

Bechtol et al.|2015l. 

ly based 
Archeo- 

For Ml 
on imaging 
logical Surv 

H, our knowledge of satellites is large 
surveys such as the Pan Andromeda 

ey (PAndAS, [McConnachie et al.||2009 

see also 

Ferguson et al.||2002||lbata et al.||2007| McConnac 

lie et al. 


20081. As shown by Brasseur et al. (20111, PAndAS is corn- 


And XXV are quiescent (Grcevich & Putman||2009 

Kaisin, 

Karachentsev & Ravindranath 

2012 Kaisin & Karachent- 

sev|2013| 

Spekkens et al.|2014 

Weisz et al.|2014a|b 

. When 


including these four lower-mass systems in our sample, the 
resulting quenched fraction remains quite high (~ 87%); fur¬ 
thermore, pushing to yet lower masses will similarly yield 
little variation in the quenched fraction, as the lowest-mass 
satellites in the Local Group are universally quenched (e.g. 
Okamoto et al.||2008| |Sand et al.||2009| |2010| |Weisz et al.| 


2014a[[Brown et al.||2014l. 


plete down to our mass limit within a survey footprint that 
covers a radial distance of roughly 150 kpc surrounding 
M31. While recent analysis of imaging data, including that 
from Pan-STARRS, has likely uncovered several new satel¬ 
lites of M31 at radial distances of > 150 kpc (And XXX- 
XXXIV, Conn et al.||2012' Martin et al.|[^13b|a| 20141, we 
exclude these systems from our study given that distances 
(relative to M31) are poorly constrained and potentially 
greater than 300 kpc. Early analysis of the stellar popu¬ 
lations in Andromeda XXXI and XXXH, however, do sug¬ 
gest passive populations devoid of young stars ( Martin et al.j 
2013bI. As likely quenched satellites, we note that inclusion 


While our census of satellites in the Local Group is likely 
incomplete at low masses due to obscuration by or confu¬ 


of And XXX-XXXIV would have little impact on the ob¬ 
served quenched fraction in the Local Group and thus on 
the results of our analysis. Moreover, while the M31 satel- 
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4 Fillingham et al. 


Galaxy Name 

Morphology 

M* (10® Mq) 

Mhi (10® Mq) 

L’hoat (kpc) 

Vlos (km s 1) 

Quenched 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Milky Way 

Sagittarius 

dSph 

21 

< 0.007 

18 

169 

Yes 

Sculptor 

dSph 

2.3 

0.22 

86 

79 

Yes 

Fornax 

dSph 

20 

0.17 

149 

-33 

Yes 

Leo I 

dSph 

5.5 

< 0.006 

258 

174 

Yes 

M31 

And I 

dSph 

3.9 

< 0.35 

58 

-82 

Yes 

And XXIII 

dSph 

1.1 

< 0.1 

127 

52 

Yes 

NGG 147 

dE/dSph 

62 

< 0.003 

142 

118 

Yes 

And II 

dSph 

7.6 

< 0.27 

184 

73 

Yes 

NGG 185 

dE/dSph 

68 

0.13 

187 

105 

Yes 

And VII 

dSph 

9.5 

< 0.36 

218 

24 

Yes 

IG 10 

dirr 

86 

50 

252 

-32 

No 

And VI 

dSph 

2.8 

< 0.015 

269 

-62 

Yes 


Table 1. Observed properties of satellite galaxies in the Local Group with measured stellar masses in the range 10® — 10® Mq and 
located within 300 kpc of either the Milky Way or M31 (separated according to host system): (1) morphological classification and (2) 
stellar mass from |McConnachie| | |2012 i ; (3) atomic gas mass me asurements from|McConnachie| l|2012[| andfCrcevich &: Putman| )|2009 I, 
plus upper limits from the literature ( Burton fc Lockman|1999| [Grcevich &: Putman||2009| [Huang et al.|[^12b[ |Giovanelli et al.|2013 l; 
(4) distance from the nearest host as given by|McConnachie| | |2012[ |; (5) line-of-sight velocity with respect to the host from |McGonnachie] 
( |2012^ , with the exception of And XXIII ( [Kirby et al.||20T^ ; (6) identification as quenched versus star-forming, according to observed 
atomic gas fraction (see Fig. |^. 


lite population may be incomplete at radial distances of 
150 — 300 kpc from M31, this should have minimal impact 
on the observed quenched fraction in the Local Group (even 
assuming that star-forming satellites are possibly biased to¬ 
wards large radial distances, as our results would imply). In 
0 we discuss this potential selection effect in more detail. 


3 SIMULATIONS 

To model the evolution of satellite galaxies in the Local 
Group, we utilize the Exploring the Local Volume In Sim¬ 
ulations (ELVIS) suite of 48 high-resolution, dissipation¬ 


less simulations of Milky Way-like halos (Garrison-Kimmel 
et al. |2014a| . The suite includes 24 isolated halos as well as 


12 mass-matched Local Group-like pairs, simulated within 
high-resolution uncontaminated volumes spanning 2-5 Mpc 
in size using a particle mass of 1.9 x 10® Mq and a Plummer- 
equivalent force softening of e = 141 physical parsecs. 
Within the high-resolution volumes, the halo catalogs are 
complete down to Mhaio > 2 x 10^ Mq, I4iax > 8 kms“^, 
Mpeak > 6 X 10^ Mq, and Vpeak > 12 kms“® — thus 
more than sufficient to track the evolution of halos host¬ 
ing Local Group dwarfs with stellar masses of > 10® Mq. 
ELVIS adopts a cosmological model based on Wilkinson Mi¬ 


crowave Anisotropy Probe 7 (Komatsu et al. 20111, with 


the following ACDM parameters: cts = 0.801, flm = 0.266, 
Ha = 0.734, Ua = 0.963, and h = 0.71. 

As hosts for the Local Group dwarfs in the stellar mass 
range of 10® — 10® Mq, we select halos in the ELVIS sim¬ 
ulations with masses of Mpeak = 5 x 10® — 6 x 10^® Mq, 


following the abundance matching prescription of jGarrison-j 
Kimmel et al. (2014aI. Abundance matching, a common 


technique for populating simulated dark matter distribu¬ 
tions with galaxies, assumes a one-to-one relation between 


a galaxy’s stellar mass and the mass of its parent dark mat¬ 
ter (sub)halo ([Behroozi et al.|2013 Moster, Naab & White 


20131. While this simple empirical approach to modeling 


galaxy formation yields great success for massive galaxies 
matching a wide range of clustering statistics as a func¬ 
tion of cosmic time (e.g. Berrier et al.|[2006 Gonroy, Wech- 


jsler &: Kravtsov||2006p , it potentially breaks down at low 
masses, where the most massive subhalos of Milky Way-like 
simulations are inconsistent with the observed dynamics of 


nearby dwarf galaxies ([Boylan-Kolchin, Bullock & Kapling- 

hat||20111|2012||Garrison-Kimmef et af.||2014b| 

Kirby et al. 

2014 Tollerud, Boylan-Kolchin & Bullock||2014 

1. 


While uncertainties remain in the potential efficacy 
or robustness of abundance matching prescriptions at low 
masses, our results are largely unaffected by changes in 
the assumed stellar mass-halo mass relation. As our fidu¬ 
cial model, we adopt the abundance matching relation of 


Garrison-Kimmel et al. {2014aI, which in our stellar mass 
range (10® —10® Mq) yields 9 subhalos as candidate satellite 
hosts within a typical ELVIS parent halo (after removing 
stripped or disrupted subhalos). When applying the shal¬ 
lower stellar mass-halo mass relation from IBehroozi et al.l 


(20131, the number of potential subhalo hosts increases ac¬ 


cordingly, but the quantitative and qualitative results with 
regard to satellite quenching remain unchanged. To account 
for potential increased scatter in the abundance matching re¬ 
lation at low masses, we also mimic scatter in the jGarrison-j 
Kimmel et al. (2014a I relation, by broadening the range of 
halo masses associated with the observed stellar mass range, 
randomly selecting 9 subhalos within 2 x 10® < Mpeak/MQ < 
2 X 10^® as satellite hosts. Again, while this biases the sub¬ 
halo population to lower masses, our results with regard to 
satellite quenching remain unaffected. 

Independent of the particular abundance matching pre- 
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scription, we account for satellites that have been tidally- 
stripped or disrupted — and thus not included in our ob¬ 
servational sample — by restricting our analysis to sub¬ 
halos for which 14iax/V"peak > 0.3, where 14iax is the cur¬ 
rent {z = 0) maximum circular velocity for the halo and 
V'peak is the value of Knax when the halo mass is at its 
maximum. While a halo’s maximum halo mass generally 
occurs prior to infall (typically at > 1.5 Rvii), V'peak is a 
good approximation for Vinfaii since a halo’s mass, in gen¬ 


eral, does not change appreciably prior to infall (Behroozi 
et al. pOli) ). By restricting our analysis to present-day halos 


with Kiiax/’l^'eak > 0.3, we eliminate all subhalos (i.e. satel¬ 
lites) that have lost greater than 90% of their mass since 
infall, thereby selecting those subhalos that are within the 
desired mass range and are highly likely to currently host 
dwarf galaxies (i.e. not disrupted or significantly stripped). 
Finally, we note that our results are not strongly dependent 
upon the adopted Knax/f^peak cut; applying a more conser¬ 
vative selection criterion (e.g. T4iax/l^peak > 0.5) does not 
significantly change our results. 

The ELVIS merger trees include 75 snapshots from 
2 = 125 to 2 = 0, providing time resolution of roughly 
200 Myr. As discussed in more detail in SQ we are in¬ 
terested in using ELVIS to explore models where satellite 
quenching occurs at a fixed physical radius from the host 
(e.g. within 50 kpc radially). The crossing time on such 
scales (~ 100 kpc), however, is typically greater than the 
median time resolution of the snapshots, such that subhalos 
would potentially pass into and back out of the quenching 
radius between two successive snapshots. To more precisely 
determine subhalo orbits within ELVIS, we map the spatial 
position of each subhalo (relative to its host) at 20 Myr in¬ 
tervals by spline interpolating the position of each subhalo 
and corresponding host halo across each of the 75 snap¬ 
shots. This interpolation also allows us to determine the 
time at which each subhalo was accreted onto its host halo 
(i.e. the infall time or tinfaii) at greater precision. Due to 
the linear spacing of the ELVIS snapshots with scale factor 
(i.e. non-linear spacing with lookback time), the difference in 
the measured infall time between our interpolated data and 
that inferred from the standard ELVIS merger trees depends 
upon the lookback time at which a subhalo is accreted. For 
systems with tinfaii < 3 Gyr, which place the strongest con¬ 
straints on the satellite quenching timescale, the difference 
in tinfaii is relatively modest, with a median offset of less 
than 0.5 Gyr for our fiducial sample of subhalos. 


4 QUENCHING MODELS 

As highlighted in ^ observations of low-mass galaxies in 
the local Universe suggest that all isolated (or “field”) sys¬ 
tems at < 10® Mq are star-forming. Thus, all quenching at 
these low masses is likely driven by environmental effects 
associated with dwarf galaxies being accreted onto more 
massive halos, such as that of the Milky Way or M31. Two 
likely quenching mechanisms operating on satellite systems 
are starvation, where future gas accretion onto the satel¬ 
lite is halted following infall into a more massive host halo, 
and ram-pressure stripping, where the cold gas in a satellite 
galaxy is violently removed via interaction with the hot gas 
halo of the host. To mimic these respective processes, we 



Lookback Time (Gyr) 


Figure 3. The cumulative distribution of infall times (tinfaii) 
as a function of lookback time for those subhalos in the ELVIS 
simulations identified as hosts of the Local Group satellites at 
10® < M*/M 0 < 10®. The solid black line corresponds to our 
default selection criteria, where subhalos are selected according 
to the abundance matching relation of |Garrison-Kimmel et al.| 
(|2014al, restricting to those subhalos with Vmax/Vpeak > 0.3 
and within the host virial radius at 2 = 0. The dashed cyan 
and magenta lines show the corresponding distribution of infall 
times when applying alternative abundance matching prescrip¬ 
tions (see ij^, while the dashed gold line traces the accretion his¬ 
tory of only those subhalos within 0.5 i?vir at 2 = 0 following the 
|Garrison-Kimmel et al.] | 2014a| l abundance matching relation. A 
quenched fraction of 90% today, as observed for the Local Group 
at these masses, corresponds to a quenching timescale (Tqueneh) 
of ~ 1.4 Gyr. 


implement two quenching models in the ELVIS simulations: 
(i) a “starvation” model in which quenching depends only 
on time since infall into the host halo and (it) a “stripping” 
model in which quenching occurs once a satellite reaches a 
fixed physical distance from the host. 


4.1 Starvation Model 


When quenching occurs via starvation, the satellite galaxy 
is deprived of any additional supply of cold gas after infall, 
such that star formation is halted following the consumption 
of its current gas reservoir — i.e. quenching occurs within 
the gas depletion timescale (ruepi = Mgas/SFR). While ob¬ 
servations suggest that atomic and/or molecular depletion 
timescales may vary with galaxy mass (Leroy et al. 200^ 


Bigiel et al.|2011||Boselli et al.|20141 |Sonierville, Popping &| 

Trager|2015| 7 we adopt a fixed quenching timescale (follow¬ 

ing infall) over the limited stellar mass range probed in our 
analysis (10® — 10® Mq). 

In the models, we define the infall time as the lookback 
time at which a halo first crosses the virial radius, Rviiiz), of 
the host. For the ~ 23% of subhalos that cross the host virial 
radius more than once, we are likely overestimating the infall 
time. However, 90% of these subhalos are accreted for a final 
time > 6 Gyr ago, thereby contributing little to no effect on 
the final inferred quenching timescale. Figure shows the 
cumulative distribution of infall times for subhalos selected 
as hosts of Local Group satellites in the simulations follow¬ 
ing various abundance matching prescriptions. Universally, 
we find that subhalos are accreted relatively uniformly over 
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Quenching Radius (Rq^ench/ Ryir) Quenching Radius (Rc^aencJRyir) 

Figure 4. Left'. The dependence of the satellite quenched fraction on the assumed quenching radius, iJquenchi fo'' our ram-pressure 
stripping model (see j )4.2[ l. The solid black and dashed cyan and magenta lines correspond to subhalo populations selected according 
to different abundance matching prescriptions (see In all cases, a quenching radius of 0.6 iJvir is required to reproduce a quenched 
fraction of 90% as observed in the Local Group. Right: The dependence of the mean quenching timescale (Tqqench) on the assumed 
quenching radius (Rquench) for the same set of abundance matching relations. For our fiducial subhalo sample, the grey shaded region 
illustrates the la scatter in the quenching timescale across the 48 ELVIS hosts. As for the starvation model, the high observed quenched 
fraction in the Local Group implies a short quenching timescale (~ 1.1 Gyr) in our ram-pressure stripping model. The typical Rvir (a-t 
z = 0) for ELVIS hosts is ~ 300 kpc. 


cosmic time, with the vast majority infalling into their host 
halo > 2 Gyr ago. 

As shown in Figure to match the observed satellite 
quenched fraction of ~ 90% in the Local Group, the starva¬ 
tion model - independent of the adopted abundance match¬ 
ing prescription - favors a quenching timescale (rquench) of 
roughly Gyr, including the error associated with the 

uncertainty in the observed Local Group quenched fraction 
(according to binomial statistics). Thus, on average, satel¬ 
lites of the Milky Way and M31 (at M* ~ 10® — 10® Mq) 
must quench within 1.4 Gyr of infalling onto their respective 
host halo. Within the individual ELVIS simulations there is 
non-negligible scatter in the derived quenching timescale, 
given the stochasticity of subhalo accretion events; assum¬ 
ing a quenched fraction of 90%, the la scatter in Tquench 
across the 12 Local Group-like host pairs in ELVIS is ap¬ 
proximately 0.7 Gyr. 

A quenched fraction of 90% in the local group could be 
dependent on completeness of the dwarf sample in this mass 
range, as discussed in the introduction. This is particularly 
true around M31, where PAndAS imaging only surveys the 
volume within a radius of ~ 150 kpc. To test the impact of 
this potential selection effect on our results, we also compute 
the distribution of infall times for those subhalos that reside 
within 0.5 i?vir today. While simulated subhalos with earlier 
infall times are slightly biased towards smaller host-centric 


molecular depletion timescales are relatively weakly depen¬ 
dent on redshift | Geach et al. |2011[[Magdis et al.|2012]|Sain-| 


tonge et al.|2013 Bauermeister et al.|2013 1, favoring shorter 
timescales at higher z. This conclusion is echoed in analy¬ 
sis of the satellite quenching timescale across a broad range 
of halo masses at intermediate redshift | McGee et al.||2011 


McGee, Bower fc Balogh|2014||Mok et al.|2014||Muzzin et al.| 


2014 Bahe & McCarthy|2015 1. Altogether, a modest varia¬ 


tion in the quenching timescale with redshift would yield no 
significant impact on our results, as our analysis is focused 
on the cumulative quenched fraction today and not on the 
detailed quenching times of individual satellites. 

4.2 Stripping Model 

To mimic quenching via ram-pressure stripping, we imple¬ 
ment a model in which satellites are quenched at a hxed 
radial distance (i?quench) from their host. In a stripping sce¬ 
nario, this radial quenching scale would be set by the volume 
over which the host’s hot halo reaches a density capable of 
stripping an infalling dwarf. To account for growth of the 
hot halo over cosmic time, we dehne our quenching radius, 
I?quench, iu terms of the host’s virial radius as a function of 
redshift (i.e. i?quench oc Rvii{z)) and require that a subhalo 
cross within the quenching radius at z < 3, allowing the 


distances (e.g. Diemand & Kuhlen 2008 Rocha, Peter & 

et al.|2005 Dekel & Birnboim|2006 1. For subhalos in ELVIS 

Bullock 2012 Wetzel, Deason & Garrison-Kimmel 20151, the 

that have crossed within the quenching radius at 2 < 3, 


correlation between radial distance from the host and infall 
time exhibits signihcant scatter, such that the quenching 
timescale derived from subhalos within 0.5 Rvii today is fully 
consistent with that based on the full sample of subhalos 
within 1 i?vir (see Fig.|^. 

In our starvation model, it is assumed that the quench¬ 
ing timescale (i.e. depletion timescale) is independent of cos¬ 
mic time. While studies at higher redshift are limited to sig¬ 
nificantly more massive systems, current work suggests that 


we classify the corresponding satellite as quenched; all other 
systems remain star-forming, independent of the time spent 
in the host halo (i.e. independent of infall time). 

Figure & shows the dependence of the satellite 
quenched fraction on the adopted quenching radius (from 
0.2 Rvir to 1 Rviv). We find that to reproduce the observed 
quenched fraction in the Local Group (90%) requires a 
quenching radius of roughly 0.6 i?vir, again with little depen¬ 
dence on the assumed abundance matching relation. Adopt- 
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ing this value for the quenching radius (-Rquench = 0.6 Rvii), 
we are able to determine the time at which each subhalo 
in the simulation (i.e. satellite) quenched following infall — 
and thus the quenching timescale (measured relative to in¬ 
fall). In Figure]^, we show the mean quenching timescale as 
a function of the chosen quenching radius (/Jquench) within 
the 48 ELVIS simulations. For large values of i?quench, each 
satellite - largely independent of orbit - spends little time 
in the host halo prior to being quenched, such that the ef¬ 
fective quenching timescale (rquench) is short. Adopting a 
quenching radius of i?quench = 0.6 i?vir, we find an average 
quenching timescale of 1.1 Gyr, in very good agreement with 


also Hirschmann et al. 20141. Interestingly, for these most 


that derived from our starvation model (see (4.11. Finally, 
a complementary analysis by Slater & Bell (20141, in which 


quenching commences at pericenter passage, finds very sim¬ 
ilar results (rquench ~ 1 — 2 Gyr, see also Weisz et al.||2015 
Wetzel, Tollerud fc Weisz||2015 1. 


5 DISCUSSION 

Given that nearly all field dwarfs are star-forming at stel- 


lar masses of < 10® M© (Grebel |l997 


Mateo 


1998 


Geha 


|et al.|2012[ ), environment-dependent processes are the dom¬ 
inant driver of quenching at low masses and must act on a 
timescale so as to avoid over- or under-producing the num¬ 
ber of star-forming satellite galaxies in the Local Group. In¬ 
dependent of the particular physical mechanism(s) at play, 
the high quenched fraction (~ 90%) observed for low-mass 
Local Group satellites thereby dictates a remarkably short 
quenching timescale, measured relative to infall. Through a 
detailed comparison of the observed Milky Way and M31 
satellites to the ELVIS suite of simulations, we show that 
this quenching timescale (rquench) is ~ 2 Gyr at stellar 
masses of 10® — 10® Mq in the Local Group. 

A very similar analysis by [Wheeler et al.| ( |2014| ), in 
which the properties of nearby dwarfs in the SDSS are 


compared to the Millennium-II simulation (Boylan-Kolchin 
jet al.|[2009p , finds that the quenching timescale must be 
much longer (> 7.5 Gyr) at slightly higher stellar masses 
(~ 10® Mq). This longer quenching timescale is dictated 
by the relatively low quenched fractions observed in nearby 
groups ( Geha et al.|2012 K While Wheeler et al. (20141 pref¬ 
erentially study higher-mass host halos than that of the 
Milky Way and M31, analysis of the Local Group is in broad 
agreement. At M* ~ 10® — 10^® M©, the satellite quenched 
fraction in the Local Group is only ~ 40% (Wheeler et al. 
|2014| [Phillips et al.|[2015| [Slater fc Bell[[2014p , which re- 
quires a quenching timescale of ~ 7.8 Gyr when compar¬ 
ing the infall times of ELVIS subhalos selected to have 
M* ~ 10®-® —10®'® Mq according to the abundance matching 


relation of Garrison-Kimmel et al. 12014a I. 


In an effort to study quenching for the most massive 


satellite systems, Wetzel et al. (20131 and De Lucia et al. 


(20121 compare group and cluster populations from the 
SDSS ( Yang et al.|2007 1 to Y-body and semi-analytic mod¬ 
els, respectively. While analysis at these masses is compli¬ 
cated by the fact that some isolated or field systems are 
quenched independent of any environment-related effects, 
these independent analyses conclude that the quenching 
timescale at ~ 10^® — 10^^ M© increases with decreasing 
satellite stellar mass, reaching as long as ~ 5 — 6 Gyr (see 


massive satellites, no evidence is found for variation in the 
quenching timescale with host halo mass ( [Wetzel et al.|20f^ 
but see also [Phillips et al.|2015 |. 

In Figure we combine our estimate of the quench¬ 
ing timescale for low-mass satellites in the Local Group 


(M* 


10'^ Mq) with the complementary constraints from 


Wheeler et al. (20141 and Wetzel et al. (20131 at higher 


stellar masses. Due to differences in the adopted definiton 
of infall time (see discussion in Wheeler et al.|2014 l, we as¬ 
sume a quenching timescale of 8 Gyr at 10® Mq. This value 
for Tquench is fully cousistent with the results from |Wheeler| 


et al. (2014 Tquench ~ 7.5 — 9 Gyr); as previously high¬ 


lighted, analysis of higher-mass subhalos in the ELVIS sim¬ 
ulations — systems corresponding to 10®'® —10®'® Mq follow¬ 


ing the Garrison-Kimmel et al. (2014a I abundance matching 
relation — yields a very similar result (rquench ~ 7.8 Gyr; 


see also De Lucia et al. 20121. Altogether, we find that 


the quenching timescale increases with decreasing mass at 
> 10® Mq, such that environmental processes are relatively 
inefficient in suppressing star formation at ~ 10® Mq. At 
lower stellar masses (< 10® Mq), however, the quenching 
timescale decreases dramatically, indicating a stark increase 
in the satellite quenching efficiency for low-mass systems. 
Given that all known satellites of the Milky Way and M31 
are quenched below stellar masses of ~ a few x 10® Mq, 
the quenching timescale is expected to remain short to even 
lower masses. For the lowest-mass systems (i.e. ultra-faint 
dwarfs) environmental quenching is likely irrelevant, as the 
smallest dark matter halos may have their star formation 
suppressed significantly by non-environmental mechanisms 


like cosmic reionization (e.g. Rees [l986 

Efstathiou 1992 

[Wyithe & Loeb|2006[ [Ofiorbe et al.|2015 

1. 


For an infalling satellite of any mass, assumed to be 
isolated from subsequent gas accretion, an upper limit on 
the longevity of any star formation activity is set by the gas 
supply within the satellite upon infall — including the satel¬ 
lite’s hot halo, from which gas may cool. Processes such as 
ram-pressure stripping, tidal effects, and harassment would 
then work to diminish this gas supply, so as to accelerate 
the quenching process. The long quenching timescales ob¬ 
served for higher-mass satellites (~ 10® — 10®® Mq) thus 
favor starvation as the dominant quenching mechanism. Ig¬ 
noring future accretion, recycling, or gas cooling, the du¬ 
ration of star formation should follow the gas depletion 
timescale (rdepi = Mgas/SFR, where Mgas is the mass of 
a satellite’s cold gas reservoir upon infall). While studies of 
nearby star-forming galaxies suggest that the surface den¬ 
sity of ongoing star formation is most closely tied to the 


dense H 2 gas (Bigiel et al. 2008 2011 Leroy et al. 20081 


observed molecular depletion timescales are much shorter 
than the inferred quenching timescales. For relatively mas¬ 
sive galaxies (M* > 10®® Mq), spanning a broad range of 
redshift, current observations point towards molecular de¬ 
pletion timescales of roughly 1 — 2 Gyr ([Bigiel et aL|[2011 


Leroy et al.|2013[|Saintonge et al.|2011[ [Tacconi et al.1|20T0 


2013[ Denzel et al.[|2010[ [2015 I. Moreover, recent work tar¬ 

geting larger samples of lower-mass galaxies conclude that 
the molecular depletion timescale decreases with decreasing 
stellar mass, such that infalling satellites at 10® Mq should, 
on average, exhaust their molecular reservoirs in less than 
1 Gyr ( Saintonge et al.|2011 Boselli et al.|2014 l. 
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Figure 5. The dependence of the satellite quenching timescale (rquench) satellite stellar mass at M* ~ 10® — 10^^ Mq. The inset 
legend lists the typical host halo mass for each data set. The magenta, gold, and cyan colored bands show the constraints for satellites 
in different mass host halos derived from analysis of galaxy groups and clusters in the SDSS tyJWetzel^^talJ_j20]^. The black square 
denotes the typical quenching timescale for slightly lower-mass satellites (~ 10® Mq) from [Wheeler et al.| ( |2014| l, with the horizontal 
error bars denoting the 25 — 75% range in stellar mass probed by that work and the vertical error bars giving the variation in 
corresponding to satellite quenched fractions of 25 — 55% (as derived from analysis of subhalo populations in ELVIS, see ijsLOur estimate 
for the low-mass satellite quenching timescale in the Local Group is given by the black circles (rquench ~ 1-4 Gyr, see 94.1|, where our 
sample of subhalos is divided into two stellar mass bins of 10® — 10^ Mq and 10^ — 10® Mq. Each point gives the quenching timescale for 
subhalos in that mass bin assuming a satellite quenched fraction of 90%, while the vertical error bars illustrate the variation in rquen^h 
corresponding to satellite quenched fractions of 80 — 100%. Our results are largely independent of mass, with increased scatter at higher 
masses due to the smaller number (and thus increased stochasticity in the infall times) of massive subhalos. The horizontal error bars 
show the quartiles of the subhalo mass distribution within each bin. At high masses, the quenching timescale increases with decreasing 
stellar mass. Below M* ~ 10® Mq, we find that satellite quenching becomes dramatically more efficient, suggesting a likely change in 
the physical mechanisms at play. 


These relatively short molecular depletion timescales 
are in direct conflict with the derived quenching timescales 
for satellites at ~ 10® — 10^® Mq, which indicate that star 
formation typically continues for > 4 Gyr after infall. Given 
that strong variation in star-formation efficiency for H 2 is 
disfavored by both observations of the Kennicutt-Schmidt 
relation (Bigiel et al. 2008 Leroy et al.||2008 Krumholz & 


McKee|2005 1 and the bimodality of rest-frame galaxy colors 

(or specific star f ormation rates, Balogh et al.|[2009 Wein-| 
mann et al.|2010| p]the extended star formation in satellite 
galaxies seemingly must be driven by a fuel supply beyond 
the molecular reservoir. While gas cooling from a satellite’s 
hot halo could replenish the fuel supply for star formation, 
there is little observational evidence for significant hot gas 
coronae around star-forming galaxies, even at the high-mass 
regime and especially when restricting to gas with a suffi¬ 
ciently short cooling time (Sun et al.|2007 [Rasmussen et al. 


2009 Grain et al.|20l6 Miller fc Bregman|2015 1. 


On the other hand, if atomic gas is included as a po¬ 
tential fuel source for star formation, there may be better 
agreement between the expected gas depletion and satellite 
quenching timescales. As previously stated, star formation 


® While a reduction in the star formation efficiency for satellites 
would yield longer depletion (and thus quenching) timescales, it 
would also lead to overproducing galaxies in the “green valley”. 


is observed to follow molecular gas density more closely than 
atomic gas density, although the weaker correlation with Hi 
surface density could be largely due to atomic gas spanning 
a broader range of densities, thus on average tracing dense 
star-forming clouds more poorly (Glover & Glark 2012a|b 


Clark et al.||2012 l. As shown in Figure [^, when including 


the atomic gas component, the inferred quenching timescales 
at high stellar masses (M* > 10® Mq) are in broad agree¬ 
ment with the measured and predicted (atomic plus molec¬ 
ular) gas depletion timescales and thus consistent with a 
picture where satellite quenching is driven by starvation at 
high masses. 

In the low-mass regime (< 10® Mq), the atomic com¬ 
ponent becomes increasingly important, with observations 
of nearby star-forming systems finding that the typical Hi 


gas fraction increases steadily towards lower mass (Skillman 


Cote & Miller 20031 Geha et al. 2006 

Leroy et al. 20081 

Schiminovich et al.|20101 |Gatinella et a . 

2010||Huang et al.| 


2012aI. More importantly, this relative growth in atomic gas 


outpaces the corresponding increase in the average specific 


star formation rate (e.g. Noeske et al. 2007 Salim et al. 2007 
Speagle et al.|20T4)), such that the gas depletion timescale is 


generally long (> 10 Gyr) for low-mass star-forming galax¬ 
ies — a result that is fully consistent with the dearth of 
quiescent field galaxies at < 10® Mq (Geha et al. 2012 


Ofiorbe et al. 20151. In comparison to the typical quench- 
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Figure 6. Left: The dependence of the satellite quenching timescale (rquench) depletion timescale (rdepi) on satellite stellar 

mass at M* > 10® Mq. The magenta, gold, and cyan colored bands show the constraints for satellites in different mass host halos derived 
from analysis of galaxy groups and clusters in the SPSS by|Wetzel et al.| l |2013| | . The black square gives the typical quenching timescale 
for lower-mass satellites 10®'^ — 10^'^ M©) from [Wheeler et al.| | |2014| |. The solid and dashed black lines show the mass dependence of 
the gas depletion timescale (H 2 + Hi) derived from combining the observed gas fractions of jBoselli et al.| l |2014| | with the star formation 
rate-stellar mass relation of [Salim et al.j | |2007[ | and [Speagle et al.[ | |2014[ | , respectively. The solid blue and solid green lines give the 
corresponding predictions from the hydrodynamical and semi-empirical models of [Dave, Oppenheimer & Finlator[ l [2011[ | and [Popping, j 
[Behroozi &: Peeples [ l (20T5| l, respectively. The inferred satellite quenching timescales at > 10^ Mq show broad agreement with the observed 
gas depletion timescales for field systems, suggesting that starvation is the main driver of satellite quenching at these masses. Right: The 
dependence of the satellite quenching timescale on satellite stellar mass at 10® < M*/Mq < 10^^, including the estimates of Tquench from 
Fig. The light grey shaded regions highlight the expected dominant quenching mechanism and projected mass dependence of rquen^h) 
while the vertical hashed grey region demarcates the critical mass scale below which a more efficient physical process (e.g. ram-pressure 
stripping) drives relatively rapid quenching in the Local Group. At high masses (> 10® Mq), satellite quenching is consistent with being 
driven by starvation; below ~ 10® Mq, however, we posit that ram-pressure stripping becomes increasingly effective, such that depletion 
(and thus quenching) timescales are significantly shorter. 


ing timescales at low masses, however, the overabundance of 
atomic gas in these systems suggests that starvation alone 
cannot drive satellite quenching at < 10® Mq . Instead, a sec¬ 
ondary process, such as ram-pressure stripping, is needed to 
remove cool gas from infalling systems, thereby decreasing 
the gas depletion timescale significantly. 


As observations of systems such as Leo I suggest | Sohn 
et al.|20T3 1, we would expect stripping to occur as a satellite 


infalls from Rvir towards pericenter, such that quenching at 
low masses would proceed according to the dynamical time 
(rdyn), which is roughly 1 — 2 Gyr for a Milky Way-like halo 
today (e.g. Stewart et al. 20091. While a detailed analysis 
of the efficacy of ram-pressure stripping (or other quenching 
mechanisms such as tidal forces) is beyond the scope of this 
work (see Fillingham et al. in prep), Figure]^ presents a 
qualitative depiction of the potential physical processes at 
play, where satellite quenching is primarily driven by starva¬ 
tion at higher masses (> 10® Mq for the Local Group) while 
ram-pressure stripping becomes increasingly effective below 
a critical mass scale of ~ 10® Mq, yielding shorter quench¬ 
ing timescales and high observed quenched fractions. In this 
qualitative picture, we expect that ram-pressure stripping 
(or some other efficient quenching mechanism) plays a crit¬ 
ical role in suppressing star formation down to the mass 
scales at which reionization inhibits gas accreti on and there- 
fore star formation (perhaps, M* ~ 10“^ Mq 


Brown et al. 


2014 Ofiorbe et al. 2015| Wheeler et al. 20151. Applying 


this qualitative picture to higher-mass host halos (e.g. rich 
groups and clusters), the critical scale at which ram-pressure 
stripping becomes effective should move to higher satellite 


masses due to the increased density of the host’s hot halo in 
concert with higher infall velocities for the satellites. 


6 SUMMARY 

Gomparing observations of the Local Group satellite popu¬ 
lation to corresponding high-resolution Y-body simulations 
from the ELVIS suite, we investigate the typical timescale 
(following infall) upon which star formation is suppressed 
in low-mass (< 10® Mq) satellite galaxies. When combin¬ 
ing our work with complementary analyses of higher-mass 
satellite populations, we present a comprehensive picture of 
satellite quenching spanning roughly five orders of magni¬ 
tude in stellar mass. Our principal results are as follows: 

• To reproduce the high fraction of quenched low-mass 
satellites in the Local Group, the typical quenching 
timescale at M* ~ 10® — 10® Mq must be relatively short 

('Tquench ^ 2 Gyr). 

• The longer quenching timescales inferred for higher-mass 
satellites are roughly consistent with the observed cold 
gas (H 2 plus Hi) depletion timescales in corresponding 
field systems; this suggests that satellite quenching is 
largely driven by starvation at M* > 10® Mq. 

• At low masses (M* < 10® Mq), the much shorter satellite 
quenching timescales are potentially set by ram-pressure 
stripping or some other process that removes the cold gas 
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reservoirs in infalling systems, such that nearly all low- 
mass satellites in the Local Group (or more massive host 
halos) are quenched. 

• If ram-pressure stripping is responsible for increasing 
the satellite quenching efficiency at low masses (M* < 
10® Mq) within the Local Group, we expect that this crit¬ 
ical quenching scale will shift to higher satellite masses in 
higher-mass host halos (e.g. rich groups and clusters) due 
to the increased density of the host’s hot halo in concert 
with higher infall velocities for the satellites. 


We have suggested here that the transition to rapid 
quenching in the smallest satellites (M* <10® M©) is po¬ 
tentially explained as a ram-pressure scale (see Figure|^and 
Fillingham et al. in prep). Above this scale, satellite quench¬ 
ing occurs on a gas depletion time, as accreted galaxies are 
starved of fresh fuel for star formation. Below this scale, 
a more rapid, drag-induced stripping becomes dominant, 
such that the gravitational restoring force of lower-mass dark 
matter halos can no longer retain cool gas against the fluid 
pressure they experience while in orbit around a larger host. 

The idea that Local Group dwarfs are quenched via ram 


pressure has been around for some time (e.g. Einasto et al. 


1974 Lin fc Faber|1983 Blitz fc Robishaw|2000 1 but the ev¬ 

idence for very rapid quenching provided here adds further 
argument in its favor, and identifies a characteristic mass 
scale where it may become dominant. Underlying this hy¬ 
pothesis is the notion that the Milky Way and M31 both 
harbor extended (~ 150 kpc) reservoirs of hot baryons with 
sufficient density to strip small galaxies of their gas. There 
is evidence that M31 does indeed host an extended circum- 
galactic medium ( Lehner, Howk fc Wakker|2015 1, and X-ray 
studies are consistent with this possibility around the Milky 


Way ( Fang, Bullock fc Boylan-Kolchin|2013 Miller Breg-| 


man 2015 Fang et al. 20151. If extended low-density coronae 


of this kind are common around ~ L* galaxies, then they 
may prove important for understanding the baryon cycle 
and the global census of baryons in the Universe. 

Of course, one shortcoming of our analysis is that we 
have restricted ourselves to satellites of M31 and the Milky 
Way. Given completeness issues, the Local Group is a rea¬ 
sonable starting point, but it remains possible that satellites 
in our vicinity are unusually quenched, perhaps owing to an 
uncommonly dense distribution of hot baryons in our vicin¬ 
ity, or atypically early infall times for the Local Group satel¬ 
lite population. Ongoing efforts to discover and characterize 
dwarf satellite systems around other massive hosts will be 
important for solidifying the existence of a sharp transition 
scale at M* < 10® Mq, where satellite quenching appears 
to become extremely efficient in the Local Group. If future 
studies reveal that the Local Group is indeed typical, then 
pinpointing the origin of this quenching scale and testing 
the ram-pressure hypothesis will be of utmost importance. 
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